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Abstract: We study weakly disordered quantum wires whose width is large compared to the Fermi wavelength. It is 
conjectured that such wires display universal metallic behaviour as long as their length is shorter than the localization 
length (which increases with the width). The random matrix theory that accounts for this behaviour - the DMPK theory - 
rests on assumptions that are in general not satisfied by realistic microscopic models. Starting from the Anderson model 
on a strip, we show that a twofold scaling limit nevertheless allows to recover rigorously the fundaments of DMPK theory, 
thus opening a way to settle some conjectures on universal metallic behaviour. 
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1 Introduction 

We discuss the heuristic DMPK (Dorokhov-Mello-Pereyra-Kumar) theory of disordered wires in 
Section [TJ following to some extent the survey by Beenakker, [2J. Then in Section |2 we introduce 
our microscopic model, a variant of the Anderson model on a strip, and we state the rigorous 
results, Propositions [4] and [5] The last proposition establishes that the main assumption of DMPK 
theory - the replacement of the microscopic model by a random matrix ensemble - is satisfied in a 
certain scaling limit. 

While we were finishing this manuscript, a related article [28] appeared on the archive. The 
results of [28 1 are very similar to ours, although the focus and motivation are different. 

1.1 Phenomenology of disordered wires 

Consider a wire of length L and cross section A, made from a disordered material and let A > be 
a measure of the strength of the disorder. We assume that electrons are injected a fixed energy E 
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and are scattered coherently in the wire. Let I be the mean free patti of the electrons in the wire 
(hence I depends on A, namely I oc A -2 ) and Xp the Fermi wavelength. We choose I as the basic 
unit of length in which L is measured and hence we define a reduced length 



L 



Physically the Fermi wavelength A p determines the number of channels N that fit into the area 
A. Here, we shall rather take the dimensionless number N as a measure of the cross section of the 
wire since, in our model, the number of channels will be naturally given. 




Figure 1: A schematic picture of the disordered wire and its various parameters. 

If N = 1, it is well-known that that the electron gets localized with localization length « I 
(hence s « 1). However, the localization length increases with N, and we can ask how the sys- 
tem behaves for s -C N , before localization sets in. This regime turns out to be experimentally 
accessible and it is called the metallic, or diffusive, regime. One of the fascinating aspects of this 
regime is the phenomenon of universal conductance fluctuations (UCF). Let g = g(s, N, A) be the 
conductance of the wire, expressed in units of the conductance quantum 2e 2 /h. It is a random 
quantity due to the disorder. Its average satisfies Ohm's law as long as s is smaller than the loc- 
alization length ~ N and it drops off exponentially once the localization length is exceeded. The 
rough behaviour for large N (see below for precise conjectures) is 

{ N/(l + s)+0(N°) s<iV (Ohm's law) 

E(g)~{ (1) 
[ cxp{— s/N} s>N (localization) 

where E(-) indicates an average over the disorder. The UCF mean that 

Var(. g ) = 2/(15/3), fors<iV (2) 

has a fixed value, independent of the material or the wire length and width. The only parameter 
that remains is the symmetry index /3 that refers to Dyson's symmetry classes. 

It is important to note that these phenomena should emerge in a large N limit only. On the 
other hand, N cannot be too large because then we enter the regime of 2-dimensional localization 
(assuming that the wire has one transverse dimension). Below, we distill some mathematical con- 
jectures from the physics literature. Keep in mind that g is a random variable that depends on s, N 
and A, and that the disorder strength A enters in the definition of the unit of length s. Moreover, 
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the conductance g can in principle depend on the way disorder is brought in and on the energy E. 
Hence, in order to state the conjectures precisely one should start from, for example, the micro- 
scopic model introduced in Section [2] (where @ = 1) and one should assume (possibly weakened 
versions of) Assumptions[T]and|2] 

Conjecture 1 (Localization lengths) For A small enough but fixed 

E(g) < const • e - "^, for alls, TV; (3) 
Conjecture 2 (Universal conductance fluctuations) 

lim lim Var( 5 ) = 2/(15/3), for all < z < 1 ; (4) 

W/oo X\0 
s/N—z fixed 

Conjecture 3 (Ohm's law and weak localization correction) 

lim lim E(g) = - + ~(l-2/P), for all < z < 1 , (5) 

JV/'oo X\0 Z 3 

s/N—z fixed 

where the constant ~ (1 — 2//?), called the weak-localization correction, is yet another manifestation 
of universality. 

At the time of writing, none of these conjectures is confirmed mathematically. In 1251 , a weakened 
version of Conjecture 1 was proven, with TV fixed and the range of admissible A's shrinking as 
TV /• oo. One should remark that Conjecture 1 differs from 2 and 3 in that the latter Conjectures 
cannot hold true without some scaling limit in which A \ i.e. I /• oo. Indeed, if we did not scale 
A, then for large s and TV one should recover the two-dimensional Anderson model, which is ex- 
pected to be localized. Conjecture 1, however, makes sense even in the regime of two-dimensional 
localizations since it states only an upper bound on the localization length, linearly growing with 
the sample width TV. 



1.2 Transfer matrices and conductance properties 

A convenient macroscopic description of the wire is through its 2TV x 2TV scattering matrix S. 
Heuristically, the S-matrix transforms incoming free waves into outgoing free waves. Let again 
TV be the number of channels and assume that the free dynamics is time-reversal symmetric; it 
follows that one can distinguish TV right-moving and TV left-moving free waves. Having this de- 
composition, the S-matrix acts as 

*($)-($) "-{*')• (6) 

where c™, c™*, c^ ut , c 1 ^ are TV-dimensional vectors representing the amplitudes of the respective 
waves. The subscript L/R indicates whether the wave is left/right moving respectively, and the 
superscript in/out distinguishes incoming and outgoing waves. The TV x TV matrices t,t' and 
r , r' are usually called transmission and reflection matrices. Current conservation is equivalent to 
the unitarity of the S'-matrix. Further, we assume that a particular basis has been chosen in 
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such that time-reversal is implemented by complex conjugation. If we assume that the scattering 
process is time-reversal invariant as well, then S equals its transpose. In short, the two symmetries 
read 

S-^S*, S = S*, (7) 

where 7 denotes complex conjugation on each matrix element. 

For our purposes, it is more convenient to encode the scattering data in the transfer matrix M, 
which transforms free waves on the left of the wire into free waves on the right. It acts as 

M ( „! ) = ( C °£]- (8) 



c. 







( 










i 



Indeed, a left-moving wave on the right-hand side of the wire, and a right-moving wave on the 
left of the wire are incoming, (they move from infinity towards the wire) whereas the two others 
are outgoing. Moreover, transfer matrices satisfy a simple multiplicative composition rule when 
two scatterers a put together in series. Finally, eqs. © and ((6) imply simple algebraic relations 
between the matrix elements of S and M , and the symmetries I0 translate into 

M*E Z M = E„ S.A/S, = M, (9) 

where we employed the 2N x 2N matrices 

S 2 = 

In other words, 

V P a 

where the matrices a, f3 are constrained to satisfy 

a*a-]?P = l, a*(3-p*a = 0. (10) 

In summary, a time-reversal invariant scattering process is described by a transfer matrix that 
belongs to the Lie group 

Q := {M G GL(2N) \ M*Yl z M = T, z and Y. X MY, X = M} . (11) 

It follows from polar decomposition of the reflection and transmission matrices (or equivalently 
of a , j3 above) that any element of Q can be written as 
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where U, V are unitary matrices and T is a diagonal matrix whose entries are the eigenvalues of 
t*t, and t is the N x N transmission matrix appearing in S, see (01. Note that the matrices U, V are 
not unambiguously determined because of the invariance of (l2)l under 

Ut+UA, Vh+AV, (13) 

where A is a diagonal matrix with the diagonal consisting of —1 and +l's. The group theoretical 
aspects of transfer matrices are discussed in detail in (24l[l8l . 

We mentioned earlier that the conductance properties of a wire with transfer matrix M can 
be read off from the eigenvalues of T. More precisely, the conductance can be defined via the 
Landauer-Buttiker formula as 

N 

g := T* = Tr T = Tr t*t (14) 

i=\ 

where T % are the eigenvalues of T. Of course, it is also possible to define the conductance through 
a physical setup, and then to prove that it is given by (fl4l >, but this will not concern us here. 

1.3 The DMPK-theory of disordered quantum wires 
1.3.1 Stochastic differential equation for transfer matrics 

The basic postulate of the DMPK theory states that the random transfer matrix A4(s) correspond- 
ing to wire length s (in units of the mean free path) satisfies an Ito stochastic differential equation 

dM{s) = dC(s)M{s), M(0) = 1 , (15) 

where C(s) is a matrix-valued Brownian motion, and d£(s) is independent of M(s) (in mathemat- 
ical terms: the process M. (s) is adapted to the filtration generated by C(s), see e.g. |20| for details). 
Equation | [T5)l is an immediate consequence of the composition rule for transfer matrices 

M{s + ds) = (1 + dC(s))M{s) 

where 1 + dC(s) is the transfer matrix for the infinitesimal piece of wire of length ds. In practice, 
one assumes that ds is large compared with the mean free path £, but small enough such for the 
transfer matrix 1 + dC(s) to be close to the identity. We immediately point out that this the latter 
requirement is not canonical. This is due to the fact that we use conventions such that, at zero 
disorder, the transfer matrix equals 10. In our mathematical treatment of a microscopic model, 
this will be assured by expressing the transfer matrix in an 'interaction picture' where the fast 
oscillations due to the nonrandom ballistic evolution are subtracted. 

Current conservation and time-reversal invariance, eqs. 10, restrict dC(s) to satisfy 

dC*{s)Y, z + H z dC(s) + d£*(s)T, z dC(s) = , 
cLC*(s)£ x + E x dC*(s) = 

as can easily be checked by the Ito calculus. In the upcoming Section II .3.21 we shall postulate 
the law of dC(s) as it is derived in the physics literature from maximal entropy considerations. 
However, for the sake of the geometrically minded reader, we remark that Huffmann [ IS J derives 
the form of the equation l(T5)l by demanding that the diffusion is generated by the Laplace-Beltrami 
operator on a certain symmetric space (note that there is no canonical Laplace-Beltrami operator 
on Q as it is a noncompact Lie group). 

4 In the terminology of |2|, we assume the 'equivalent channel' assumption rather than the 'isotropy' assumption 
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1.3.2 The "Maximum entropy" assumption 

With l fl5)l at hand, the major assumption of the DMPK theory is the 'isotropy', or 'maximal en- 
tropy' (MEA) assumption which states that dC(s) is drawn from an ensemble of maximal entropy, 
constrained on the total scattering, Tr(l — T(dC(s)), due to d£(s). We do not go into this derivation 
and we merely state its conclusion, namely that C(s) is the following matrix- valued process, 



where 



bij(s) = 



co = ( ! (s) b(s) 

' { b(s) a(s) 



l/V2N-{BK(s)+iB^(s)) 

i/y/N-B^s) 

-Oji(s) 

1/V2N- (£«(*)+ 
bji(s) 



(16) 



1 < i < j < N 
i =3 

otherwise 

1 < i < j < N 
otherwise 



(17) 
(18) 



and BR(s) , B^ (s) , B 1 ^ (s) , Bfjj(s) are independent standard Brownian motions for 1 < i < j < N. 
We note that by construction 

a = —a* , and b* = b , 



and we have the following invariance property: For any unitary matrix W, 



WaW* 
WbW* 



(19) 
(20) 



where ~ denotes that the random variables have the same law. In particular, this invariance prop- 
erty means that we need not specify the basis in I|17lll8|l . Hence, ia is GUE-distributed, as one 
could also have noticed from the explicit expression (fl7b . As regards dC(s), our choice for o and b 
implies that 

WdC(s)W* - dC(s) (21) 



where 



W 



W _0_ 

w 



— I for any unitary W. 



(22) 



This expresses the intuition that the disorder in the wire is not 'biased' towards any specific chan- 
nel basis. In conclusion, the SDE l fl5l l, together with the isotropy assumption encoded in the law 
defined by j2"Tll2"2"l l, defines a stochastic flow in the group Q (MEA-flow). 

Assume that we aim to calculate the transmission eigenvalues Ti(s) associated to the transfer 
matrix M.(s). This can be done by Ito calculus, since the eigenvalues Ti(s) are Ito processes that 
are functions of M.[s). We postpone the explicit result of this calculation to the next section. First, 
we explore the main consequence of the symmetry j2"TT l of dC(s). 

We write T(s) = T(A4(s)) for the diagonal matrix of transmission eigenvalues associated to 
M(s) by the decompositon ((12). Likewise, we write U(s),V(s), remembering that they are not 
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uniquely determined by M(s). We also abbreviate 



U{s) := 
V(s) := 
Af(s) := 



U(s) _0 

U(s) 

V{s) 

V(s) 



V^TR v^Hi) - 1 



Let W be as in l l22t with W Haar-distributed in the unitary group. Then, formally, 

T(M{s + ds)) = T((l + d£{s))U(s)Af(s)V(s)) 

~ T(W(1 + d£(s))W*U(s)N{s)V(s)) 

~ T((l + d£{s))W*N(s)V(s)) 

~ T((l + d£{s))W*Af(s)) . 

To obtain the equivalence on the second line, we used (|2"T)l . and for those on the third and fourth 
line (which are actually equalities) we observe that W* (s)U(s) ~ W* (s) since the unitary operators 
W in W(s) are Haar-distributed, and that the diagonal matrix T(s) is not affected by the unitaries 
in the decomposition lj!2|l. Importantly, this computation shows that T(s + ds) depends on A4(s) 
only via T(s) (because it only depends on Af(s)), whereas in general, the eigenvalues at s + ds ought 
to depend on both the eigenvalues and eigenvectors at s. 

M(s) — s- M(s + ds) 

MEA-flow 

I I (23) 

T{s) — > T(s + ds) 

D M P K — equation 

This phenomenon, represented by the commuting diagram (23) , is known in statistical mechanics 
as the autonomous evolution of a (set of) macroscopic degree(s) of freedom. Of course, it is rather 
well known in random matrix theory as well, as it appears in Dyson's Brownian motion, see e.g. 
(8lH7). Another conclusion is that, for the sake of calculating T(s + ds) from T(s) = T(A4(s)), we 
could as well have assumed that the unitary U(s) in U(s) is Haar-distributed. We can also turn the 
previous argument around, to obtain the following statement 

If, for all s, the matrix U(s) that appears in the polar decomposition dl2l of M (s) is Haar distributed in 
the unitary group and independent ofT(M(s)), then T(M(s + ds)) can be calculated form T(M(s)) as if 
dC(s) were distributed as in ([16}. 



This observation suggests a heuristic explanation as to why the DMPK equation can still be 
a good approximation when dC(s) is not distributed as in dl6t . By using the 'concentration of 
measure' property on the unitary group, see e.g. [17J, we understand that the set of U(s) such that 
the dependence of T(A4(s + ds)) on T(M(s)) deviates significantly from that predicted by the 
DMPK equation, has an exponentially shrinking size w.r.t. the Haar measure, as N /* oo. 
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1.3.3 The DMPK equation 

In the previous section, we already anticipated the fact that the evolution of the transmission 
eigenvalues 71- (s) is autonomous. We now give the precise evolution equation for them. 

Since T k (s) = T k (M.(s)) is a function of M(s), the law of its process can be calculated by Ito's 
fomula: 

BT 1 B^T 

^(-) = ESA^W^(-)+2 E dM jM^ )dM ^ a)dM - {a) 

The derivates g M k . , gM ° , t g M .. can be computed by spectral perturbation theory under the as- 
sumption that the diagonal matrix T(s) is non-degenerate, i.e. all transmission eigenvalues are 
simple. 

The formal result is the DMPK equation (discovered by II71 H91 ') 

dT k (s) = v k (T(s))ds + D k (T(s))dB k (s) (24) 

where B k (s),k — 1,...,N are independent standard Brownian motions and the drifj| v k (T(s)) 
and diffusion constant D k (T(s)) are given by (we keep the parameter 8 to make the comparison 
with the literature simpler, in our case (3 = 1) 

2T k I T k + Tj — 2T k T,j 



v k = -T k + x 1 - T k + 8/2 V ■ 



/3N + 2-/3 



We refer to [2] for a more detailed account of this straightforward, but somewhat lengthy, calcu- 
lation. We also note that the term 'DMPK equation' usually refers to the Fokker-Planck equation 
associated with the SDE l|24|l . The Brownian motions B k (s) in the equation | |24|| originate from 
linear combinations of the Brownian motions in $Y7\ . We call the result j24l l formal since we have 
not established that the equation admits a solution for all times, i.e. we have not excluded that two 
or more eigenvalues can collide, see [17J on a possible strategy to do this. 

The basic picture concerning |(24| is that, as s grows, all T k (s) are driven to by the term — T k in 
the drift (note that < T k < 1). However, due to the repulsion in the term containing (T k — Tj)^ 1 , 
the Tfc's keep a distance of order 1/7V. Finally, for small values of T k , Tj, the repulsion vanishes 
and the T k 's pile up at (in the localized regime, they are all very close to 0). The most interesting 
question is now whether the equation (|24) l reproduces the phenomena discussed in Section fTTTI In 
other words, 

Does the random variable g satisfy the conjectures listed in Section \Ll\ if one defines g — J2i T% with Ti 
the solutions of the equation j24|l ? 



5 We follow the literature on the subject in calling and D k the drift and the diffusion. One could however argue 
that this is a confusing convention, since absence of drift, v k = does not render the process time-reversal invariant, cfr. 
difference between Ito and Stratonovich conventions. 
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On the heuristic level, the answer is clearly "yes" for all 3 Conjectures, see e.g. [2J for an ex- 
pansion scheme for large N. On the rigorous level however, the situation is not so clear (at least 
to us). For Conjecture 1, the techniques of [23), developed for a harder problem, are sufficiently 
strong to settle the question, see e.g. (22). The perspectives seem to be best for (3 = 2 where the 
DMPK equation can be solved exactly [3J in terms of toroidal functions. It is plausible that in this 
case, one can prove all Conjectures listed, but we do not know any reference where this is actually 
done. For (3 = 1, 4, the DMPK equation was solved in [4, 18, 30 ], but in those cases, the solution is 
less explicit. 

2 The microscopic model: Anderson model an a strip 
2.1 The Hamiltonian for a disordered wire 

Our system is an infinitely extended wire of width N , which we model by the Hilbert space Z 2 (Z x 
{1,...,N}) = l 2 (Z) <g> C N . We consider the motion of a single electron through the wire. It is 
governed by a Hamiltonian consisting of a ballistic part, and a disorder term, H\ = Hbai + XV, 
where A governs the strength of the disorder. The ballistic Hamiltonian of the system is given by 
a longitudinal Hp acting only on l 2 (Z), and a transversal H±, acting only on C N ; 

flhal = <8> 1 + 1 ® flx, (25) 
and we will slightly abuse the notation by writing 

#11 = J2j| ® 1, H X = 1®H±. 
We choose H» to be the standard lattice Laplacian on l 2 (Z), with kernel 

Hp[x,x') = £ x , x /_i + 5 x , x '+i, x,x' 6Z. 

The disorder is modeled by a random on-site potential, which we assume to be located only on 
the sites with longitudinal coordinate between 1 and L. Hence, L is the length of the disordered 
part of the wire and 

x=l,...L 
z=l,...N 

with l( Xt g) the one-dimensional projector on the site (x, z) € Z x {1, . . . ,N}, and v XfZ are i.i.d. 
random variables with mean and variance 1. We also assume that the third, fourth and fifth 
moments of v x ,z exist, i.e. 

E(v)=0, EO 2 ) = 1, E(v 3 ),E(v 4 ),E(v 5 ) < oo. (26) 

The operator H± g B(C N ) can for some purposes by chosen largely at will (though we will 
always impose Assumptions [T] and |2]below) but for our main result we will require that H± is the 
transverse Laplacian with periodic boundary conditions and with a magnetic field that breaks the 
chiral symmetry (to be introduced precisely later). The eigenvectors of H±, denoted by i/j^, \i = 
1, . . . , N fix a basis in C N , 

flxW. = -Ej-MVv.. i*=l,...,N, (27) 
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I.i.d. random variables v. 




L 

Figure 2: The microscopic model 



with corresponding energies E±(h)- The index \i is hence a free channel-index. 
After diagonalizing the free transverse dynamics, we define 

E\\ := O* (E -H 1 _)0 = diag(£|| (1), . . . , E\\ (N)) , 

where the eigenvalue -En (/i) = E — E±(fi) represents the longitudinal energy in the /ith channel. 
The obvious condition for a channel to be conducting is that Eu (fi) lies in the spectrum of Hn : 

Assumption 1 (Elliptic channels). Recall that the spectrum of Hn is the interval [—2,2]. We demand 

E\\ (h) G int(specif|| ) for all (i = 1, . . . ,N . 

This assumption does for example exclude a Hu describing a barrier that would stop the waves, 
even in the absence of disorder. 

For reasons that will become clear later (see Section I3I21 , we need to make sure that the system 
has no accidental symmetries. A first requirement would be that the operator H± is nondegener- 
ate, but we demand a stronger condition, namely: 

Assumption 2 (No degenerate level spacings). Let Hj^ be nondegenerate, i.e. Eu (//) ^ En (v) for 
H 7^ v. Let 9^ = 9(E\\ (h)) be defined by 

2cos^ = ^||(/i), < < 7T . (28) 

Let m G {1, . . . N}for i = 1, 2, 3, 4. Then the equality 

4 

X^w =0 > where q t G {-1,1} 

i=l 

can only hold (up to a permutation of the indices {1,2,3, 4}) if 

Hi = ^2 and us — Hi an d qiQ2 = ©94 = -1 • 

As Eu (h) is the energy available for longitudinal propagation in the /ith-channel, 9^ corres- 
ponds the absolute value of the longitudinal momentum in each channel. Eq. { |2"8l is indeed noth- 
ing else than the dispersion relation for Hu . 
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2.2 Transfer matrix 



The setup described above suggests to study the (random) scattering matrix S = S(X, L, N) of the 
wire, since the motion is ballistic outside a compact region. 

Abstractly, the scattering operator is obtained by comparing the full dynamics to the free one 
generated by the ballistic Hamiltonian d25t , 

S = a- lim Uo(-t-)U(t--t + )U (t+), (29) 

t±— >±oo 

where U(t) = cxp(-Hxt) and U (t) = exp(—H ba it). The question of asymptotic completeness, 
namely the existence of the limits and the unitarity of the scattering operator, is readily solved by 
noting that the perturbation XV is of finite rank and invoquing the Kato-Rosenblum theorem [21 [. 
Finally in the fiber decomposition H = J %EdE, the scattering operator is described by a family 
of 2N x 2N scattering matrices Se- 

If Assumption [TJ is satisfied, the fiber spaces T-Le are isomorphic to C 2N = C^f ® C^, which 
physically corresponds to a left and a right moving sector. Once the incoming energy is fixed 
and the index E dropped, this definition coincides with the more physical one, eq. 10, based on 
transmitted and reflected plane waves, see e.g. [1 [. In that representation, current conservation is 
equivalent to the unitarity of the scattering matrix, which is ensured by asymptotic completeness. 
If, in addition, the matrix S is written down in the eigenbasis {ip^ ® , ffi Vv}^=i defined in (j27j, 
time-reversal invariance indeed corresponds to S* = S with 7 the complex conjugation on each 
matrix entry. 

As noted in Section 1, we shall work in the following with the equivalent transfer matrix M = 
M(X,L,N). In order to construct it explicitly, we first introduce the family of matrices T x , x — 
1 , . . . , L, which can be constructed without explicit reference to scattering, and this is the way this 
setup is usually presented in the mathematical literature. Define the 2N x 2N matrix 

-^)=(0'«W( E - J f- (30) 

where H xx := l x Hl x = H± + XV X and V x = Yl x =i v x,z^( x ,z)- The key property of T X A is that, if 
^ = (^^ez, £ is a (not necessarily normalisable) solution of the eigenvalue equation 
i7 A * = E^l, then 

( \ rp\ f * x \ 

\ m x ) ±x \ J ■ 

Of course, this property is preserved under multiplication, i.e. 

( ^ ) = T\L) [%\\ T\L) := T L X-i ■ 

To describe our results in the most natural way, it is convenient to express the transfer matrix 
in the basis of free in- and out- states, i.e. in left- and right-moving waves in the different channels 
H, as was done in Sections ll.2l and [1.3l Therefore, we introduce the transformation T mapping the 
transfer matrix T x in the position representation to the more physical transfer matrix 

A4 A = T- 1 2*T, x = 1,...,L. (31) 
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For N = 1, this change of basis amounts to 



7)_r(J), (-T)-C 

(the vectors on the left should be thought of as incoming and outgoing 'plane waves') and for 
N > 1, we have N copies of this transformation, with e = Eh (/z) now channel-dependent, see j45t 
for an explicit expression. Note that although it represents a mere basis transformation, T is not 
unitary since the bases are not orthogonal. 

In that representation, the transfer matrix M A is diagonal in the absence of disorder, namely 

M° = Ml = ( exp ^ «) . ° 

\ ex P( _1 ^(^||)) 

and corresponds indeed to its heuristic definition. In the following lemma, we summarize the 
precise relation between the microscopic model and the matrices M A . This is just a mathematical 
confirmation of the discussion in Section [T] 

Lemma 3. The transfer matrices M A , defined through (f30|) and $3l\ belong to the group Q, see ftl) . They 
are related to the scattering matrix, as defined in | |29] |, by the relations given in Section \L2\ 

Note further that the transfer matrices constructed in this section represent the influence of one 
lattice site. The transfer matrices for longer parts of the wire are introduced below. 



2.3 Results 

2.3.1 Scaling limit 

Naively, we interpret the DMPK theory as suggesting that the transfer matrix for the wire with 
rescaled length s, 

M A (LA- 2 S J):=Mf A _ 2sJ ...M 1 A , 

(where [a\ stands for the largest integer not greater than a) converges to the solution of the dif- 
ferential equation | fl5l l as A \ 0. This cannot be correct as such. First of all, A/ A ([A~ 2 sJ) contains 
some rapidly oscillating terms (as functions of A~ 2 s) whose phase cannot be well-defined as a 
function of s. Recalling that the transfer matrix in the absence of disorder (V = 0) is of the form 

ex P (iLA- 2 s J0(£||)) \ 

exp(-iLA- 2 s J )) J ' 

we introduce a first modification to our convergence conjecture. We multiply by the inverse of this 
free scatting matrix to subtract the fast oscillations, i.e. we consider 

A x {[\- 2 s\) := (Af°(LA- 2 sJ)) _1 M A (LA- 2 sJ) (32) 

which is 1 at V = 0. The new conjecture is that A x ( [A~ 2 sJ ) converges to the solution of the DMPK 
equation. This weakened form is still not correct without further assumptions although we believe 
that it is 'essentially correct' (as we shall explain below). Actually, A A ([A _2 sJ) converges to the 
solution of a SDE that differs from the DMPK equation because it has a lower symmetry. 



M°(LA- 2 S j) = 
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Proposition 4. If Assumptions\T\and\2\hold, then A X ([X 2 s\) converges in distribution^ 

A x {[\- 2 s\) ^> A(s) /oralis. 

distr 

The process A(-) is the strong solution of the SDE 

dA(s) = dZ(s)A(s), A{0) = 1 



(33) 



(34) 



where 



Z(s) 



a' (a) b'(s) 
b>(s) a'(a) 



The matrix-valued Brownian motions a'(s) , b'(s) are defined as 



-a' 



1 < i < j < N 

i = 3 
otherwise 



KM 



j 1/V2N- (b^{s)+iBI u {s) 
(the various standard Brownian motions appearing here are similar to those of (T7\\TM ) with 



1 < i < j < N 
otherwise 



M 2 



NE\V^\ 2 
4 • | sin 9^ sin 9 V 



(35) 



(36) 



(37) 



where E|V^„| 2 denote the covariances of the N x N matrix V = 0*VO of the random potentials in the 
channel basis and the momenta 9^ were defined in Assumption^ 



We note immediately that the SDE j34|l differs from the corresponding SDE (TT5b in DMPK 
theory by of the factors that do depend on the channels. For the same reason, the definitions 
05l l and d3"6t are truly basis-dependent (p. v index the channels). This remains true even if the 
Hamiltonian H± is chosen in the most isotropic way as outlined below. 

Let H± on C N be given by the kernel 



H ± (z, z') := h± (c'-'S^,^ + e-*S ZtZ , +1 ) , z, z> e Z/NZ 



(38) 



(with < 7 < 2ir/N and h± > 0) which corresponds to a cylindrical wire permeated by a magnetic 
field 7 along its symmetry axis. Note that we cannot choose 7 = 2irk/N, k € Z, without violating 
Assumption |2 In those cases, H± corresponds to the pure Laplacian on a discrete torus, which 
has a chiral symmetry and hence twice degenerate energy levels. The choice d38b also forbids the 
injection energy E to lie at the band center. Indeed, since the spectrum of H± is symmetric around 
zero for any 7, the condition E = implies that the longitudinal energies also come in opposite 
pairs and therefore the momenta in pairs (9^, tt — 9^), which also violates the Assumption. 



6 As visible from the proofs, the convergence actually holds on pathspace equipped with the Skorohod topology. This 
obviously implies convergence of single-time distributions 
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Proposition 5 (Convergence to the DMPK process). Let H±_ be of the form d38l> with 7 and hj_ such 
that Assumptions\l\and\2\hold. Let A(s) = A(s)' 1 ' h± be the process defined by Proposition^ Choose a 
sequence h±(n) such that for any n, Assumptions\l\and\2\are satisfied and h±(n) \ as n /• 00. Then 

^( s )T,fcx(«) "Z^° M{s/c), (39) 

distr 

where M{s) is the solution of the DMPK-SDE l fl5l l and c — c(E) is a constant that only depends on the 
energy E. 

Note that the limit in | |39l becomes independent of 7. Finally, by combining Propositions |4] and 
[5j one gets 

Urn lim ,{A\[\- 2 S \)) =M(s/c), 

h ± (n)\0 X\0 

namely the convergence of the microscopic Hamiltonian model to the DMPK random matrix 
model. The order of limits cannot be exchanged, since the bare limit h±(n) \ eliminates the 
transverse motion and hence all interaction between the channels. In fact, if h± = 0, then the 
model consists of N copies of a one-dimensional chain and hence its localization length is inde- 
pendent of N. The rescaling constant c = c(E) originates from the fact that our length unit A~ 2 is 
not precisely equal to the mean free path; see the footnote on Section ITTT1 

3 Discussion 

3.1 Related work and outlook 

3.1.1 Scaling limits in the Anderson model 

The existence of the metallic regime in the Anderson model in dimension d > 2 is a major open 
problem of mathematical physics. An important step ahead was taken with the establishment of 
kinetic E5l [13 ] and diffusive ITTTI scaling limits, where time and space are rescaled by appropri- 
ate powers of the coupling strengh A, as the latter goes to 0. In those scaling limits, the unitary 
evolution goes over to the linear Boltzmann equation or the diffusion equation. In spirit, our 
result is comparable to this approach, although the quasi one-dimensional metallic regime that 
we study is distinct from the one in d > 2. In particular, the quasi ID metallic regime exists in 
the scaling limit only (see the remarks following the Conjectures in Section fTTTb . We have a kin- 
etic space rescaling x H> A~ 2 x, in which the Hamiltonian model with disorder reduces to a fully 
stochastic model, just as the unitary evolution reduces to the fully stochastic Boltzmann equa- 
tion in the works mentioned above. To make our results complete, we should establish that the 
stochastic models obtained in Propostions|4]and|5]do possess the properties of the metallic regime. 
This seems straightforward for the scaling limit ((39b , since the eigenvalue process associated with 
this SDE is the DMPK equation, which has been studied extensively in the physics literature, see 
Section [1] In fact, the random matrix ensemble that solves the SDE ( Tl5t is often mentioned au 
par with the more familiar Dyson classes of random matrices, see e.g. [5|. To establish that the 
SDE obtained in Proposition [4] also has the universal properties of the metallic regime seems a 
more challenging task. In some sense, this is comparable to proving that e.g. Hermitian Wigner 
matrices with fully broken time-reversal invariance share the universality properties of GUE, as 
was recently established in |[T0l[T2l . 
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3.1.2 Alternative models for the DMPK theory 



Several models have been constructed to provide an alternative justification for the DMPK theory. 
Dorokhov [6] has constructed a Hamiltonian model of N chains coupled only by random hopping 
acting strictly transversally. For each pair of chains (/x, v), the hopping potential U )lv {x) at point 
x is a Gaussian variable of zero mean and independent of all other pairs. The model is therefore 
invariant under the orthogonal group by construction. In fact, it is very similar to Wegner's N- 
orbitals model [29J. Another model that reproduces the DMPK theory is the (ID) supersymmetric 
nonlinear a model, [9). 



3.1.3 Lyapunov exponents 

In 1251 , Schulz-Baldes calculated the lowest Lyapunov exponents for the Anderson model on a 
strip, at small disorder. As already remarked, this yields a partial confirmation of Conjecture 1 
in Section ITTT1 The drawback of this approach is that, by their very nature, Lyapunov exponents 
do not reveal much information on the metallic regime of disordered wires, although one could 
argue that the conjectured equidistancy of the Lyapunov spectrum is a sign of universality. In a 
recent work, the authors of [24 1 conjecture (and confirm numerically) that a certain property, the 
Random Phase Property (RPP), holds for the Anderson model on a strip. They remark that the RPP 
is weaker than the maximal entropy assumptions (see Section [l.3.21 made in the DMPK theory, yet 
it allows to estimate the Lyapunov exponents and hence the localization lengths. It is not clear to 
us (but it does not sound unreasonable) that the RPP in fact also implies universal behaviour in 
the metallic regime. Further, it seems plausible that the transfer matrix ensembles that we obtain 
in Propostion|4]satisfy the RPP. 



3.1.4 Hyperbolic Brownian motion 

The DMPK equation |24} for N = 1 was solved explicitly in [16]. It can be related to Brownian 
motion in the hyperbolic plane. A very similar description can be given from the point of view of 
the Brownian Carousel |27|. As already mentioned, the DMPK equation for any TV can be related 
to Brownian motion on symmetric spaces, see 118. 1301 . 



3.2 Idea of the scaling limit 

We now explain heuristically why the propositions hold true. Recall the definition of A x (x),x e N 
in j32"l l, i.e. A x (x) := [M°]~ x M x (x). For convenience, we define G such that 

mo ° 

-9(E\\ 



M u = exp{iG}, hence G -- 
and we drop the A-dependence on ^4 A (a;). Then, A(x) satisfies the stochastic difference equation 



A{x + 1) - A(x) = [e-^+^GTT^J^c' 10 - ll A(x) . (40) 



■-xz x 



We remark that the matrices Z x+ i are independent of A(x) and E(Z x+ i) = 0, hence A{x) is a 
discrete- time martingale. The Z x are however oscillating in x due to the factors c lxG . We can write 



e 



ixG R x e ixG , \R X := e _iG TT*T _1 - 1 
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where R x are i.i.d. matrix-valued random variables. From the explicit formula for T^, we can now 
check that the expression is indeed linear in A. Note also that the 2N x 2N matrix Z x contains only 
N independent random variables, namely the v x ,z, z = 1, . . . , N. 

To gain some insight, we first replace Z x by R x . That is, we consider the process D{x) defined 

by 

D(x + 1) - D(x) = XR x+1 D(x), D(Q) = 1. (41) 
First we remark that, with appropriate conditions on the random variable v x . z , we have 

A V ^> TZ{s), for alls, (42) 

^ — 4 distr 
x=0 

where the matrix-valued Brownian motion IZ(s) is defined by replacing each random variable v XyZ 
in the definition of R x by the Brownian motion B z (s). Here z is the index ranging from 1 to N. 
Of course, l|42) l is nothing more than a multidimensional version of the convergence of random 
walk to Brownian motion. We emphasize that the 2N x 2N matrix entries of IZ(x) are a linear 
combination of N independent Brownian motions only. It is not hard to believe that the discrete- 
time process (|4T1 l converges to the solution of the corresponding SDE 

D{[\- 2 s\) T>(s), where dV(s) = dK(s)V(s), X>(0) = 1. 

distr 

as is easily proven by standard martingale theory, see e.g. [14J and Section|U 

We now look for the analogue of the convergence (@2j upon replacing R x by Z x . Let us first 
observe it, for example, at the level of the second moment. Indeed, 



E 



E Z A * E 

x=l \ x' = \ 



kl- 



e-isCGi-Gfj+Gk-GOj 

x=l 



(Rx)„ (Rx) 



(43) 



where the indices i,j,k,l refer to a basis in which the matrix G is diagonal and the numbers 
Gi,Gj, . . . are the eigenvalues. This is hence the basis in which the free transfer matrix is diag- 
onal, and the base vectors can be indexed by the double index (/i, q) where /j, G {1 , . . . , N} indexes 
the channels and q £ {+1, —1} is +1 for the left moving and —1 for the right moving sectors. The 

eigenvalues are Gj = q9^ if j = q). Since E (R x )^ {Rx) kl is independent of x, ll43l converges 



to 



Thus, one can check that 



kl 



LA" 2 *] 

\ Y' Z x ^> Z(s), for alls, (44) 

^ — ' distr 



where Z(s) is a matrix valued BM characterized by the covariances 

E{z ij ( s )z kl (s)) = s G ^ Gj ^ Gk E(n tJ ( s )n kl ( s )) 

E (Z* (s)Z,,(s)) = 8 Gi _ GhGl -G k E (n* 3 { S )1l kl { S )) 
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It is in this place that we need to use Assumption[2] Indeed, Assumption[2]states that the eigenval- 
ues of the matrix G do not have any type of degeneracy hence a condition like (Gi — Gj+Gk—Gi) — 
is only satisfied when i = k, j = I or i = j, k = I . In other words, the fast oscillations originating 
from the free part of the transfer matrix kill most of the correlations between matrix elements of 
Z(s). As a result, the number of independent random variables in Z{s) is 0{N 2 ) whereas in IZ(s) 
it was O(N). In |28| , this phenomenon is called 'noise explosion'. 

Next, we calculate the nonvanishing covariances of Z(s) arising from our microscopic model. 
The matrix R x is given explicitly by 



XR X = e~ iG T(0* ® 1 2 ) 
= — Ae~ iG Y 



E — H x 
In 

Y — i 



-Ijv 




\ (O® 1 2 )T- 1 -M° 



v x 





where V x = V* was defined in Proposition [4] Hence 



Vx o 




T" 



Recalling the definition of the map T, we can write it explicitly as a matrix with diagonal N x N 
blocks. 



(45) 



M£||))- 1/2 


M£||))- 1/2 


(ip^r^expH^i)) 


(ip(Sn))- 1 / 2 exp(i0(£||)) , 



where p(-Eii) = 2sm9(E\\). By direct calculation, we can now check that the covariance of Z{s) 
agrees with that of the PJ-1S of (|37t , i.e. 



n(R x h. 



We now explain heuristically how this gives rise to Proposition [4] Recall that A(-) satisfies the 
stochastic difference equation 



A(x + 1) - A(x) = \Z x+1 A(x) 



(46) 



Since Z x has zero mean and unit variance, it takes a a time of order A~ 2 for A(-) to change appre- 
ciably. However, on this timescale, the oscillations in Z. are not longer visible and one can hence 
replace \Z\x-* a \ ~ dZ(s), i.e., by a time average on the fast time-scale, in |@6) . This suggests that 
the scaling limit A s of A(LA~ 2 sJ ), if it exists, must be a solution of dA(s) = dZ(s)A(s), and hence 
Proposition |4]holds. 



4 Proof of the scaling limits 

Proposition |4] is essentially a textbook result in homogenization theory, where one starts from an 
evolution equation with a fast degree of freedom that evolves independently of the slow degrees. 
In our case, the fast degrees of freedom are the oscillating phases in the random variables Z x . For 
a treatment of these matters, we refer the reader to fl5l , where one discusses a setup identical to 
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our case (see Chapter 6 and 7 of H5]). In fact, the only difference is that our model is defined on 
the lattice instead of in the continuum (and hence we can have independent random potentials 
instead of a rapidly decaying dependence in the x variable). One can mimick the proof of [15 [ 
apart from the fact that we have to establish tightness of processes with sample paths in D R d [0, oo[ 
instead of C R d [0, oof. 

For the sake of explicitness, we outline a proof based on an expansion. We first show that the 
second, third and fourth moments of A x converge to those of A(s) in a scaling limit. Then, in 
Section l4~2l we invoke a standard result to eventually prove in Section |4~3l that the convergence of 
moments is sufficient to conclude the proof. 



4.1 Convergence of moments 

We will establish the convergence of fourth moments of .A A (|_A _2 sJ ) by a brutal, but completely 
standard estimate. We use the assumption that the third and fourth moments E(v 3 ), E(v 4 ) are 
finite. 

Lemma 6. Let P e ,e — a,b,c,d be operators on C 2N with \ \P e \\ = 1 and let s e € [0, r], e = a, b, c, dfor 
some finite r. 

Let (A x (x))# stand for either (A x (x))* or A x (x). Then 

[TrP a (^ A (LA- 2 Sa J))#P fc (A A (LA- 2 Sb J))#P c (A A (LA- 2 Sc J))#P (i ( J 4 A (LA- 2 Sd J))#] (47) 



E[TrP a (A(s a ))*P b (A(s b ))*P c (A(s c ))*P d (A(s d )) 



#1 



<C(X,t,N) (48) 



where C(A, t, N) \ as A \ 0. 



Proof. To avoid too many constants, we will treat the case where s e = t, P e = 1 and # are chosen 
such that we take the adjoint of the first and third A. The general case is completely analogous. 
By iterating the difference equation j40)l , we obtain 

m = e ( Tr n (i + az*j n (i + \z Xb ) n (i + \z* Xc ) n (i + \z* Xd ) \ m 

\ X a — 1 Xb = l X c — 1 Xd — 1 / 

where we have abbreviated t\ := [A _2 rJ and the Z Xc are ordered with decreasing indices and the 
Z* have increasing indices. Evaluating the products yields 

= 1 

+ E A " E ETi ' (( z k ■ ■ ■ z k-) ( z *:> ■ ■ ■ z *i) { z k ■ ■ ■ z k-) ( z *7* ■ ■ ■ z ^) ^ 



n>0 P n ^~ 



In (|50|l , 1 < x\ < . . . , x™" < t\ are sites (elements of N), for e = a,b,c,d and some numbers 
m a + nib + nic + md = n. The sum J2 Pn runs over all such sets of sites, such that the E(-) is not 
zero, i.e. such that for each x\ that appears in the set, there is at least one partner x\ 1 , — x\ with 
e! 7^ e. This allows us to divide the sites that appear in the sets into three classes: those that appear 
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twice, three or four times. Let n%, 713, 77,4 be the number of such sites. Then 2n 2 + 3n 3 + 4774 = n. 
The number of P n corresponding to (n.2, 713, 774) is estimated as 



#{P n ■ (n2,n 3 ,n 4 ) is fixed} < 



(t a C 2 4 )" 2 (r A C|)"3 (r A C 4 4 )"* 

■ r?.o' T1.a\ 



n 2 \ n 3 \ n 4 

where C 4 3 4 are binomial coefficients. For each choice of P n , 

\r(P n )\ < {NE(v 2 )) ri2 ■ {NE(v 3 )) n3 • (ATE^ 4 ))™* . 
This gives us an a-priori bound on the sum in (50) : 

$>"5>(P„)|< J2 A 2 ^+3«3+4n 4 J2 |r(P„)| 

n Pn n2,«3,«4>0 Pn-( n 2,n 3 ,ni) fixed 

(ClE(v 2 )NT X ) n2 {ClE(v 3 )NT X ) ns {CfE{v A )NT X ) ni 



n2,n3,n4>0 



ni\ 713! 724! 



< 



E 



(const ■ Nt) U2 (const • AiVr)" 3 (const • A^iVr)' 



n 2 ! n 3 ! n 4 - 

no ,n3,n4 >0 

<e const Arr , for any A 

Hence, by dominated convergence, we can interchange the sum over 77 and the limit A \ to 
obtain 

= 1 + y lim A™ V r(P n ) 

n>0 P„ 

One can also see immediately that all terms with 773 7^ or 774 7^ vanish as A \ 0. Therefore, 
the sum over P„ reduces to a sum over pairings, each involving 71 = 2v matrices Z™ , and the 
expectation value is a simple product of v covariances. Each such pairing can be described by first 
specifying a 'pattern', namely from which group a, b, c, d the matrices are drawn in each pair, and 
then the particular point x associated to each pair. Given a pattern, we can use the same oscillatory 
argument that led to j44)l to handle the limit of the sum over all sets of points satisfying that pattern, 
thereby replacing sums over Z" L " by stochastic integrals over Brownian motions dZ(s" lc ) (recall 
that the Kronecker S's are included in the definition of Z(s)). Hence, 

limA 2y Vr(P 2l ,)= V 

E[Tr [ dZ*{sl)...dZ*(s^)- [ dZ{s™*)...dZ{sl) 

■ [ dZ*(sl)...dZ*(s™°). I dZ( s ^)...dZ( s 1 d )) . 

J0<sl<...<S™ c <T J0<S 1 d <...< S J d <T ) 

As an illustrative example of the discussion above, we perform the computation explicitly in the 
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case m a = nib — 2, m c = rrid = 0. Then 



A 2 r(P 2+2+0+0 ) = A 2 J2 H^iKlKl- Z *l Z *0) 

1,2 1,2 
x a i^f, 

\<X 1 <X 2 <T X 



A\0 



E ( Tr f dZ{s 1 )* dZ{s 2 Y dZ{s 2 )dZ{s 1 ) ) (51) 



= E I Tr / / dZislYdZislr •// dZ(s 2 )^(4) (52) 

y Jo Jo jo jo y 

where dBll follows from the convergence of the second moments, eq. ll44t , and the Gaussian prop- 
erty of the Brownian increments, and l |52l by their independence for unequal times. We note 

r— F=rn 

that the contraction pattern Z* t Z* 2 Z x iZ x x does not appear because of the ordering of the matrix 

a a ^ 

product. Moreover, the iterated Ito-integrals are well-defined since J Q a dZ(s a ) is an Ito-process 
which is adapted to the filtration associated with the Brownian motion Z(s 2 a ). 

To conclude, we observe that the summation over n yields the stochastic Neumann series as- 
sociated to the SDE j34l in each of the four brackets, and therefore ( |48i . □ 



In fact, one could repeat this proof for all moments if one assumes that all moments of v x>z 
exist. 

4.2 An abstract result on convergence to diffusions 

Let X n (t),n e N, be a sequence of IR^ -valued discrete time martingales, defined via the discrete 
time stochastic difference equation 

X n (t + A 2 ) - X n {t) = A„ • e{t)X n (t), t e A 2 N , 

where A„ \ as n /*■ oo, and = (£,ij(t))t j=i are random variables with mean zero and 

covariance 

E(er i (*)&(*')) = kt'C^ k M *J e a 2 n . 

We assume that the functions C"j k ; (•) are uniformly bounded and 

f dsf(s)C? Jtktl (s) n ^ Oijou [ T dsf(s), for any JeC^). (53) 
Jo Jo 

The jump process X n is viewed as a process in continuous time whose sample paths belong to 
D M d[Q, oo[, the space of right-continuous functions with left limits, equipped with the Skorohod 
topology. The set of jump times tj is thus A^N. In this context, we shall denote by X n (tj—) the 
left limit of the sample path of X n at t j. 

The following lemma is a version of Theorem 4.1 in Chapter 7 of [14J, simplified such as to 
match the setup above. 
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Lemma 7 (Convergence to diffusions). Let X be a R d -valued diffusion process X defined as the unique 
strong solution of the SDE 

dXi(s) = ^2 <7ij x j( s )dBj(s), X(0) = c e R d , 



with respect to the filtration generated by the d-dimensional Brownian motion B(s). Assume that 
• The processes 



W%{s) :=Y,C? M {s)X%{s)X?{s) 



k.i 



satisfy 



sup 

0<t<r 



J* ds I WTjis) - |] u^iX^Xf is) 



distr 



0; 



• The jumps of X n become small in the sense that 

E(sup tje[0>T] |Xf(t J ) - X?(tj-)\ 2 ) 
E(sup tje[0 , T] |H^(^)-WS(O-)|) 

where the sup ranges over all jump time tj. 
Then, 





0, n /* oo 



distr 



x?{t)xm)- / dswpjs) 



(54) 



(55) 
(56) 



on the space D & d [0, oo [. 
Note that 



are martingales with respect to the filtration generated by X n (t), as required by Theorem 4.1 in 
Chapter 7 of d. 



4.3 Proof of Proposition ID 

To prove Proposition HJ we use Lemma [7]with X n (s), resp. X{s), the Revalued processes con- 
taining all real and imaginary parts of the matrix elements of A x (\\~ 2 s\ ), resp. A(s), (hence 
d = 2 (27V) 2 ) and with n labeling a sequence X n such that A„ \ as n / co. The functions 
Cfj k ;(•) are the oscillating covariances of the random variables Z x , see Section 13^21 

Lemma [6] yields that all 4th moments of X n converge to those of X as n /• oo, uniformly on 
compacts, i.e. 

sup [E(XZ( Sl )X? 2 (s 2 )X? 3 (s 3 )XZ(s 4 )) -E (X tl ( S i)X t2 ( S2 )X t3 (s^XM)] -^0, (57) 

{si}e[T,T'] 
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and we note that the same argument can be repeated to provide the convergence of 2nd and 3rd 
moments. 

We check condition (54l l of Lemma [7| First, by the definition of Y n , and using the uniform 
bound on the fourth moments of X n , we observe that for any t ,t' € [0, r] and for alH, j = 1 , . . . , d, 

E\Y£(t')-Y£(t)\ 2 <(t-t'fC, 
and C < oo can be chosen independent of n. The Chebyshev inequality then reads 

P(|W)-1$(*)| 2 >&) < l(t-t'fO- (58) 

To show that this yields the convergence $54$ , let us assume that the random variable there does 
not vanish in distribution. That implies in particular that there is a sequence of times t m 6 [0, r], 
and some constants c, d > such that; 

P ( | Yg [t m ) | > c) > c', for all m . (59) 

Let t be an accumulation point of the sequence (t m ). Then 1 1591 implies that either P(|F i ™(t m ) — 
Y if(t)\ > c/2) > cf/2, which is excluded by ggj, or P(|^ (t) | > c/2) > c'/2. But the latter can also 
be excluded by observing that, as m — > oo, 

Ed^r^i 2 ) — > o- 

This follows by the convergence of fourth moments, the continuity of sample paths of X and the 
weak convergence l l53l . Hence, condition j54|l is proven. 

We now turn to condition (55). Recalling that X n (tj) - X n (t.j-) = X n ■ £, n (tj)X n (tj) for any 
jump time t j, we have for any a > and b > 0, 

P{sup tje[0 , T] \X?(tj) - X?(tj-)\ > b) 

<\- 2 -supP(\\ n -^(tj)X?(tj)\>b) 
tj 

<\- 2 - S up±E\\ n -Q(tj)X?(tj)\ a 

< i • A- 2 + a -sup(E|^.(t J )r«) 1/o (E|X;(t J )| Q!a ') 1/a ' ! 1/a+l/a' = 1 

where the first equality follows from the fact that there are A~ 2 jumps in the interval [0,r], the 
second and third are Chebyshev and Cauchy-Schwartz inequalities, respectively. If a, a' are chosen 
so that aa' < 4 and E|£y (t j) \ aa < oo, the last expression is bounded by 

Let now b > 1. To conclude, we compute 

E(sup tje[0iT] \X?(tj) - Xl l (tj-)\ 2 ) 

OO 

< +F(\H n \ g [A^,l]) + ^E(|^"| 2 || J ff n | G [6 fe ,6 fc+1 ]) .P(|£r| G [fe fe ,6 fe+1 ]) 

oo 

< + 0(A^ ) + 0(A- 2+a ) ■ b ' 2[k+1) ■ b ~ ka , 

k=0 
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where we have abbreviated H n = sup tje m T i \Xf{tj)—Xf{tj—)\. Condition j55l l follows by noting 
that the choice a > 2 ensures both the finiteness of the series and the convergence of the expecta- 
tion value to zero. Finally, Condition | |56|| follows by a similar argument. This concludes the proof 
of Propositions] □ 
Note that a > 2 implies a > 2 and therefore aa > 4, whence the need to control higher 
moments of the potential v, see j26l l. 



4.4 Proof of Proposition [5] 

To prove Proposition [5j we note that, for H± as in (f38b - the basis of eigenfunctions f/v consists 
simply of the exponential functions on Z/NZ, independently of 7. Hence, 

E I %v 1 2 = 1 /N, independently of fi, v 

Since hx \ 0, also E±Xp) \ for all fi and hence E\\ (p) -» E and M -> where is the solution 
of 2 cos 8 = E. This means that 



v 4(1 -(cosfl) 2 ) 4(1 -(E/2) 2 ) ' 

which fixes the scale factor c in the RHS of (|39l l. Since both the processes .4. and M. are solutions 
of SDE's with smooth and uniformly bounded coefficients, convergence of the solutions follows 
from convergence of the coefficients, see e.g. Chap. 11 in l26l . □ 
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